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We propose and quantitatively develop two schemes to quickly and accurately generate a stable initial con- 
figuration of neutral atoms in optical microtraps by extraction from the Mott insulator state in optical lattices. 
We show that thousands of atoms may be extracted and stored in the ground states of optical microtrap arrays 
with one atom per trap in one operational process demonstrating massive scalability. The failure probability 
during extraction in the first scheme can be made sufficiently small (~ 1CT 4 ) to initialize a large scale quantum 
register with high fidelity. A complementary faster scheme with more extracted atoms but lower fidelity is also 
developed. 



PACS numbers: 03.67 .Lx 

Introduction: Recently quantum computation using neutral 
atoms has attracted much attention because neutral atoms are 
well isolated from the environment and may offer a route to 
scalable quantum computation. Different neutral atom quan- 
tum computing proposals |H 0. [H differ in trapping tech- 
niques. We focus on arrays of microscopic traps (micro- 
traps) that can be implemented using independent focused 
laser beams. The typical size of a microtrap in current ex- 
periments is about 2pm because of the diffraction limit of a 
laser beam as well as experimentally technical issues. Unlike 
an optical lattice, all microtraps can be moved independently 
in position space, therefore they have advantages in single 
atom addressability and controlled interactions between pairs 
of atoms. In recent years, many schemes have been proposed 
for universal quantum computation in microtrap systems JH]. 

These schemes are based on the assumption that thousands 
of neutral atoms can be prepared in the ground states of optical 
microtrap arrays with one atom per trap, i.e. they assume the 
ability to initialize a neutral atom quantum register. Although 
impressive experimental progress has been made in trapping 
single atoms [4], such an assumption has not been fulfilled be- 
cause the trapping processes are random and not determinis- 
tic. Furthermore, the trapped atoms are not in the trap ground 
states. These difficulties have prevented neutral atom quantum 
computation architectures from accomplishing the kind of im- 
pressive experimental progress recently achieved in ion trap 
quantum computation J5|]. To overcome these difficulties one 
proposal seeks to extract single atoms from a Bose-Einstein 
condensate (BEC) by moving an optical dipole trap out of the 
condensate |g]. 

In this paper, we propose two concrete schemes to quickly 
extract thousands of neutral atoms from the Mott insulator 
(MI) state in optical lattices [ 7] and prepare them in the ground 
states of optical microtrap arrays with precisely one atom per 
trap. In the MI regime, atom number fluctuations at each site 
are suppressed and a MI state may be obtained with proper ex- 
perimental parameters, as demonstrated in recent experiments 
yl 0, [Hi- Indeed, a perfect MI state with one atom per lattice 
site may be implemented through purification processes @]. 
In MI states, atoms are isolated from each other and occupy 
the ground states of each site making it difficult to realize a 
low density MI in long lattice spacing (d) laser geometries 
ifloll where the MI energy scales (~ d~ 2 ) fall below charac- 



teristic temperatures. On the other hand, the size of a micro- 
trap (~ 2pm) is much larger than the typical lattice spacing 
(~ OApm), which make it difficult to extract atoms to micro- 
traps directly from a MI state. Therefore the main challenge, 
and also the goal of this paper, is to transfer single atoms from 
optical lattices with a short lattice period (SPOL) to spatially 
separate optical microtraps. 

We show in the first scheme that such transfer process 
may be accomplished using hyperfine state dependent opti- 
cal lattices with a long lattice period (LPOL), microwave ra- 
diation, and resonant "removing" lasers. A LPOL is cre- 
ated by intersecting two laser beams at a certain angle, as 
demonstrated in a recent experiment 1 1 Oil - The LPOL in- 
duces position-dependent energy shifts of the hyperfine states 
of atoms, which, when combined with microwave radiation 
and resonant "removing" lasers, expel many atoms out of the 
optical lattice and form a patterned loaded optical lattice with 
one atom per n lattice sites (« > 3 is an integer). The re- 
maining atoms in the lattices are well separated, and can be 
adiabatically transferred to the ground states of optical micro- 
trap arrays with one atom per trap. With our scheme, thou- 
sands of atoms may be extracted and the failure probability in 
extraction can be kept below ~ 10 4 . In the second scheme, 
the transfer process is accomplished using state-dependent fo- 
cused lasers without expelling atoms out of the optical lattice, 
therefore more atoms can be extracted but the failure probabil- 
ity is much higher. These complementary schemes, one with 
very high fidelity (but relatively lower speed) and the other 
with very fast speed (but relatively lower fidelity), take ad- 
vantage of the robust localization inherent in the MI to isolate 
atoms. The whole operational process is simple and within 
currently accessible experimental technology. 

The first extraction scheme includes four operational steps 
and their time sequences are schematically plotted in Fig. la. 
In the following, we explain each step and study various er- 
rors that could lead to failures in extraction. For simplicity, 
we focus on a one dimensional geometry but we emphasize 
that our technique can be straightforwardly applied to two di- 
mensional arrays. 

Step (I) Initial Mott state: Consider a pure 87 Rb BEC pre- 
pared in the hyperfine ground state |0) = \F = I, trip = -1) 
and confined in a quasi-one dimensional (x direction) har- 
monic magnetic trap. Along the transverse direction, the 
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Figure 1: (color online). Schematic plot of the single atom extrac- 
tion, (a) Time sequence for laser depths at four extraction steps. 
Solid line (SPOL), dotted line (LPOL), dashed dotted line (remov- 
ing laser), dashed line (microtrap laser), (b) Geometry of additional 
lasers defining the LPOL. (c) Site specific energies in the optical lat- 
tice. Solid and dashed lines correpsond to the optical potentials for 
atoms at states |0> and |1) respectively. A labels target atoms for ex- 
traction and B labels other atoms. 



atomic dynamics are frozen out by high frequency optical 
traps it 1 111 . An optical lattice along the x direction with wave- 
length A s = 850«m is ramped up adiabatically (~ 200ms) to 
a large potential depth of V s = 50Er such that the BEC is 
converted into a MI state with roughly one atom per lattice 
site, using properly chosen trapping parameters and number of 
atoms 10]. Here Er = h 2 /2mA 2 denotes the recoil energy. This 
MI state may have defects [8] and clever purification schemes 
(|9|] may yield a nearly perfect MI state with exact one atom 
per lattice site. Such a perfect MI state is the reservoir for 
single atom extraction and provides the starting point of our 
scheme. 

Step (II) Selective Depopulation: Here atoms at specific, 
unwanted sites are transferred to another hyperfine state. Two 
<x + -polarized laser beams intersecting at an angle 8 (Fig. lb), 
which drive the 5S — > 5P transition and has a detuning 
A2 = -2n x 3608G//z to the 5 2 P3/2 state (corresponding to 
a wavelength A/ = 787. 6nm), are adiabatically ramped up. 
These two beams form a LPOL along the x direction with the 
lattice period r\\ - A// [2 sin (6/2)]. Here we require 77/ to be n 
times that of the short wavelength optical lattice period A s /2, 
that is, 6 = 2 arcsin [Ai/nA s ]. The LPOL induces energy shifts 
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for two finestructure ground states |±) = 
|5S : j = l/2,mj = ±1/2), where I(r) is the intensity of 
the laser, T q (q = 1,2) is the decay rate for states 5 2 Pi/2 
and 5 2 Pj/2- cj q (A q ) is the frequency (detuning) for the 
transition from 5S to 5 2 Pi/2 and 5 2 Pi,/2- c+i = 0, c + 2 = 1, 
c_i = - V2/3, c_2 = VI /3 are transition coefficients. Two 
hyperfine states |0> and |1) = \F = 2, m = -2) can be written 
as |0> = 1/4 1-> + 3/4 |1> = |~) using the respective 
Clebsch-Gordan coefficients. Therefore the energy shifts 
for states |0> and |1> are 6E = 6E-/4 + 36E+/4, and 



5E\ = 5E-, which give the total shift of the hyperfine splitting 
between |0> and |1>, SE = 6E X - 6E = - where 

a q = 3nc 2 r q I/2a>q. Following a similar procedure, we find 
the spontaneous emission rates for atoms at |0) and |1) are 
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Notice that the wavelength A/ = 787.6nm (i.e. the detuning 
A2 = -27rx3608G//z) is optimized to obtain the maximal ratio 
77 = 6E/y between the shift SE and the rate y = max \jq, y{\. 

The energy shifts induced by the LPOL are spatially depen- 
dent for atoms confined in the SPOL. In Fig. lc, we plot the 
combined lattice potentials for atoms at |0) and 1 1) with n — 3. 
We see that the LPOL induces two different shifts of the hy- 
perfine splitting, depending on the positions of atoms (A or 
B), and the difference 6 = 6E (A) — 6E (B) can be adjusted 
and chosen to be 6 — 52Er. Applying the adiabatic condition, 
we estimate the ramp up time to be 44/js that corresponds to a 
10~ 4 probability for excitation to higher bands. 

We then apply a microwave n pulse to flip the quantum 
states of atoms at position B (Fig. lc) from |0) to |1). The 
microwave is resonant with the hyperfine splitting of atoms at 
B, but has a detuning 6 for atoms at A. A time-dependent mi- 
crowave pulse with frequency Q(f) = Qoexp^-w 2 ,? 2 ) (— tf < 
t < tf) is used to perform the n pulse. For a set of mi- 
crowave parameters wo = 6/4 — \3Er/U, tf = 5/ojq, and 

Q = it I 
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23Er/H, the pulse flips the 

quantum state of atom B from |0) to |1) in 2tf = 38.6/js, while 
the error to flip atom A to state |1) is found to be 5.9 x 10~ 6 by 
numerically integrating the Rabi equation lfl2Tl that describes 
the coupling between |0) and 1 1 ) through the microwave pulse. 
The LPOL are adiabatically turned off after the microwave 
pulse. During the whole process time r, the probability for 
spontaneous scattering of one photon from each atom is esti- 
mated to be P = C ydt » 1 x 10~ 4 . 

Step (III) Remove Non-target Atoms: In this step, atoms at 
B are removed from the trap by applying a cr^-polarized "re- 
moving" laser that drives a resonant cycling transition |1) — > 
|2) = |5 2 P3/2 : F = 3,nif = -3J. Scattering photons from the 
laser push (or heat) non-target atoms B at |1) out of the trap 
without affecting target atoms A at state |0) because of the 
large hyperfine splitting (v m 2k x 6.8G//z) between the two 
states. To push (or heat) atoms out of a trap with depth Uq, 
the number of spontaneous emission photons needs to be at 
least n p — Uo/2Er |4]. For instance, 25 photons are needed 
for each non-target atom to remove it from an optical lattice 
with depth Uq - 5QEr. The dynamics of the photon scatter- 
ing process are described by the optical Bloch equation lfl2tl . 
from which we can numerically calculate the number of scat- 
tering photons n p for both target and non-target atoms. We 
find the number of scattering photons can reach 25 in a short 
period ~ 1/xs for atoms B, but it is only ~ 10~ 5 for atoms 
A. Therefore the impact of the resonant laser on the target 
atoms A can be neglected. Note that hot non-target atoms 
tunnel more easily in the optical lattice, which enhances the 
collision probability between target and non-target atoms at 



different lattice sites. However, because of the high la 
depths (~ 50E R ), the tunneling rate is quite low (~ 0.01 
on average for hot non-target atoms, which yields a long 
neling time (~ 100ms). The short lifetime (~ lps) of 
non-target atoms makes the collision probability very si 
(~ Ips/lOOms = 10~ 5 ). The effect of interatomic collis 
on the target atoms can therefore be neglected. 

Step (IV) Transfer to Microtraps: In this last step, the 
maining target atoms A at the optical lattice sites are a 
batically distributed to the ground state of optical micrc 
arrays with one atom per trap. The optical microtraps 
be focused far red-detuned lasers or high frequency b 
detuned optical traps using Hermite-Gaussian TEMoi rr 
beams 11X111 - Each optical microtrap is focused near one 
get atom and contains only one atom because of the 1 
spacing between atoms in the lattice. In the distribution 
cess, the effective potentials the target atoms experience 
be approximated as harmonic potentials with the trapping fre- 
quency vy(f) = J{4Vf (t) /w 2 + 2V L (t) k 2 ) Jrn, where V f (t) 
and Vl (?) are the potential depths of the microtraps and the 
optical lattice respectively, w is the beam waist of the micro- 
trap lasers. Vf (?) and Vl (0 may be varied simultaneously to 
adjust m (t) from its initial value m (0) = (2Er/K) VVl (0) JEr 
to the final expected m (to). If m (to) = tu (0), the final mi- 
crotrap potential depth Vf (to) = V L (0)k 2 w 2 /2 * \366E R 
(~ lOApK) for initial parameters V s (0) = 50£r, and w = lpm. 

When a deeper or shallower microtrap potential depth is 
needed, the trapping frequency m(t) should be adjusted so 
that the adiabatic condition is satisfied 
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where |</>^, and \<p e ) are the ground and the excited state wave- 
function, AE g is the energy gap between two states, and £ is 
the adiabaticity parameter. Because of the parity of the wave- 
functions, the lowest possible excitation is to the second ex- 
cited state, which gives AE g = 2hzn, and {<j> e \ |2 = %/ V2. 
The adiabatic condition, Eq. ([31, yields a time-dependent trap- 
ping frequency m(t) = m(0) / (l + 4 V2£cr (0) tj, where + 
correspond to a deeper and shallower final trapping frequency 
respectively. 

Under the adiabatic approximation, the quantum states of 
the target atoms can be expanded using the time-dependent 
basis <p(f) = c g (t)(j>g(v7(t))+c e (t)(p e (:vj(t)), where <p g (w (i)), 
and (p e (xu (t)) are the adiabatic ground and second ex- 
cited states of the Hamiltonian H (m (?)) with the associated 
eigenenergies E g - hm(t) /2 and E e = 5hm(t) /2. Inserting 
this expansion into the the Schrodinger equation for a single 
target atom A yields a coupled equation for the coefficients 
c g (t) and c e (t): 



dt \ c e 
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where x(t) = -ih^ <0cl ||l fe > = %AE g . This equation can 
be solved analytically to give the occupation probability at the 
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Figure 2: The adiabatic process for the distribution of the target 
atoms from the optical lattice to arrays of optical dipole micro- 
traps with higher trapping frequencies. Adiabaticity parameter is 
£ = 0.005. (a) Trapping frequency with respect to time, (b) The 
probability P e (t) = |c e | 2 for atoms to be at the second excited state of 
the effective harmonic trap with respect to time. 



second excited states 

P e (t) = \£ sin 2 (in ( 1 + 4 V2£?zr (0) f) /4 V2^) . (5) 

In Fig. 2, we plot the trapping frequency m(t) and excita- 
tion probability P e (t) with respect to time for a deeper final 
trap. We see P e (t) oscillates with time, but has the maxi- 
mum 4£ 2 . For adiabaticity parameter £ = 0.005, the maxi- 
mal excitation probability is 10~ 4 . The total distribution time 
is T = (1 -m(Q) lm(T)) I(A^tj(Q)). For instance, T ~ 
9Aps is needed to transfer target atoms from the optical lattice 
to microtraps with trapping frequency m(T) = 4?zr(0) (cor- 
responding to a 1.66 mK trap potential) with 10~ 4 excitation 
probability. We note that this transfer time is much shorter 
than the characteristic hopping time (~ 5s for V s = 50Er) of 
atoms within the depleted optical lattice. 

Combining all four operational steps, we find that thou- 
sands of atoms can be extracted from the optical lattice to mi- 
crotraps in less than 300/is. For instance, in a one dimensional 
optical lattice with 300 usable atoms for single atom extrac- 
tion, we can extract 100 atoms (one per three lattice sites). In 
a two dimensional trap, the same process can extract 1/9 of 
total atoms (~ 10 4 ) simultaneously, which offers a large scale 
up in the initialization stage of a quantum computation using 
microtrap arrays. 

The failure probability of extraction does not decrease with 
an increasing n. That is because the microtrap potential is al- 
ready very weak at a ~ 1. 3pm displacement (the position of 
the neighboring atoms in a n = 3 superlattice) from the trap 
center for a typical ~ 2pm microtrap. The potential has a neg- 
ligible effect on the process of transferring neighboring atoms 
to other microtraps. A large n lattice does not contribute to the 
largest error source of the scheme: the heating due to sponta- 
neously scattered photons from atoms in the state-dependent 
lattice in step (II). 

A Speedup Scheme: In the above scheme, the qubit sup- 
ply time for quantum computation is limited by the period 
for preparing Bose-Einstein condensates, which is typically 
on the order of minutes. Notice that in step (III) most atoms 
(2/3 in ID or 8/9 in 2D lattices) are lost from the lattice by 
applying the resonant "removing" laser. In the following, we 
propose a scheme for single atom extraction without remov- 
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egy balances these two sources for the degradation of initial 
qubit fidelity. Assume that the optical potentials along the y 
direction for target atoms are 
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Figure 3: (color online) (a) Schematic plot for the single atom extrac- 
tion from paralell one dimensional lattices. The distance d = 5fim. 
Filled and opened circles correspond to target and non-target atoms 
respectively, (b) Optical potentials for target atoms along the y direc- 
tion. The length and energy units are a c and h 2 /2mcrl respectively. 
a(t) = 0.2<r c (solid line); a(t) = 0.8cr c (dashed line); a(t) = 1.5<x t . 
(dotted line). 



ing non-target atoms, which can then be recycled for further 
extraction processes. In this scheme, 87% of the atoms can 
be extracted within several superfluid-insulator transition cy- 
cles (~ Is), considerably speeding up the rate for supplying 
fresh qubits for a neutral atom quantum computer. However, 
the scheme induces much larger spontaneous emission prob- 
ability of photons and excitation probability to high motional 
states for atoms, which significantly degrades the fidelity of 
the initial qubit, therefore further cooling is needed to obtain 
qubits with high fidelity. 

A schematic for the speed up scheme is plotted in Fig. 3a. 
Consider a series of parallel well separated (d = 5 pm) ID 
optical lattices (along the x direction) in the xy plane. In the 
speedup scheme, steps (I)-(II) are still preformed to place all 
target atoms for extraction into state |0) and non-target atoms 
to |1). We then ramp up a focused laser with the detuning 
A = -2n x 780GHz from the 5S -> 6 2 P 3/2 transition so 
that the focused lasers only induce red-detuned traps for tar- 
get atoms at state |0), but do not affect non-target atoms at 
state |1) lfl3ll . The focused lasers adiabatically move along y 
direction to take target atoms out of the optical lattices, with- 
out affecting non-target atoms at state |1) (Fig. 3 a). The atoms 
inside the focused lasers are then adiabatically transferred to 
far-detuned dipole traps to suppress spontaneous emission of 
photons. These dipole traps are moved along the x direction 
to a cooling bath to improve the fidelity of the initial qubits. 
At the same time, the Mott insulator states in the optical lat- 
tices are melted by adiabatically ramping down the depth of 
the optical lattice. The trap parameters are adjusted so that an- 
other Mott state with one atom per lattice site can be obtained 
as the lattice depths are adiabatically ramped up. Repeating 
the above steps, we can extract another 1/3 of the remaining 
atoms. After 5 such cycles, about 1 - (2/3) 5 ~ 87% atoms 
can be extracted. Finally, the remaining atoms inside the opti- 
cal lattices are discarded and a new BEC must be produced to 
continue the process. 

To extract the atoms adiabatically from optical lattices 
avoiding excitations to higher bands of the focused laser traps, 
high potential depths are needed, which leads to high sponta- 
neous photon scattering probabilities. Therefore a good strat- 



V(y) = -V c exp(-2{/ 2 /cr 2 ) 

-V f exp(-2(y-a(t)) 2 /afj, 



(6) 



where V c and V/ are the potential depths for the ID con- 
finement and the focused laser, respectively. Gaussian beam 
approximations with waists cr c and cry have been used for a 
rough estimate, a (t) is the position of the focused laser cen- 
ter and its rate of change should satisfy the adiabatic con- 
dition h \da (t) /dt\ = £A£ 2 / \{<p e \ dH/da \p g ) , where the en- 
ergy gap AE g and transition matrix element {(f> e \ dH/da \<f>g)\ 
can be evaluated for different potentials V (a (f)) by solv- 
ing the single particle Schrodinger equation. The poten- 
tial minimum position y m \ n for the target atom is deter- 
mined through dV/dy = 0. The wavefunction of the tar- 
get atom is expanded in a harmonic oscillator basis *T„ (y) = 
(x/ti) 114 exp(-x(y- y min ) 2 /2}H„(^y) / V2"n! around the 
potential minimum y m \ a , where the oscillation frequency xr = 
(d 2 V/dy 2> j . /m and H„ ( yfxyj is the Hermite polynomial. 
The Hamiltonian for the target atoms evaluated in this basis 
for different a (t) gives the matrix representation H„ m (a (f)) = 
0?nl H I*!*,,,) (n, m < 10 is enough for our calculation). Diag- 
onalization of the matrix yields the eigenenergies and eigen- 
functions, which determine AE g and \{<f> e \ dH/da \<pg^ for dif- 
ferent a (t). The total moving time can be estimated using 

naf 

T= h\{(f> e \dH/da\(/> g )\/tAE 2 da (7) 
Jo 1 ' 

where aj is the final position of the focused laser. 

In Fig. 3b, we plot the optical potential V (y) for a set of 
parameters cr c ~ 0.93/um, crj ~ 0.46//m, V c = 200h 2 /mcr 2 , 
Vf = 280ft 2 /ma 2 and three different a (f). Applying the above 
procedure with these parameters, we estimate the extraction 
time ~ 5ms and the excitation probability to high bands ~ 
7 x 10~ 3 , with the spontaneous scattering probability ~ 10~ 2 . 
We see that the fidelity of the atoms is not as high as that for 
the original scheme. Therefore further cooling is needed to 
significantly improve the fidelity of initial qubits. 

Conclusion: We propose two schemes for extracting thou- 
sands of atoms simultaneously from the Mott insulator state 
in optical lattices to optical microtrap arrays with one atom 
per trap. The extracted atoms stay at the ground states of the 
microtraps. In the first scheme, about 11% of atoms from 
a BEC can be extracted and the failure probability is ~ 10~ 4 
with properly chosen experimental parameters. In the speedup 
scheme, about 87% of atoms from a BEC can be extracted, but 
the failure probability is much higher ~ 10~ 2 . We provide a 
detailed quantitative analysis validating the feasibility of our 
proposed schemes for neutral atom quantum computation. 

We thank S.L. Rolston for valuable discussion. This work 
is supported by ARO-DTO, ARO-LPS, LPS-NSA. 
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